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Supersymmetric gauge theories on the lattice 

I. Montvay a 

a Deutsches Elektronen Synchrotron DESY, 
Notkestr. 85, D-22603 Hamburg, Germany 

The perspectives of numerical simulations in supersymmetric quantum field theories with vector-like gauge 
symmetries are discussed. A numerical simulation algorithm for SU(2) gauge theory with gluinos is studied and 
the first results on the glueball-gluinoball spectrum are presented. 
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1. INTRODUCTION 

Supersymmetric gauge theories may have fas- 
cinating non-perturbative properties. Previous 
attempts to understand them theoretically are 
mainly based on the instanton calculus (for a re- 
view and references see Recently, the N = 2 
extended gauge theories received particular inter- 
est because Seiberg and Witten were able to ob- 
tain exact results on their low energy effective 
action [||. Besides the theoretical interest, there 
are also some interesting open questions concern- 
ing phenomenology. Although the most probable 
range of supersymmetry breaking mass param- 
eters, as for instance the gluino mass, are at or 
above the 100 GeV range, the possibility of a light 
gluino with a mass near 1 GeV is not yet com- 
pletely excluded || . In such a situation the extra 
hadronic states ( "R-hadrons" ) could be hidden in 
the complicated "zoo" of normal hadronic states 
with masses above 1 GeV [@. 

In a lattice formulation, most probably, one has 
to give up exact supersymmetry. This is certainly 
not surprising, because supersymmetry is an ex- 
tension of the Poincare symmetry, which is bro- 
ken by the lattice. For instance, the supersym- 
metry generators Qi a ^Q^ with internal indices 

i,j and Weyl-spinor indices a,/? (i,j = 1,..,N; 
a, (3 = 1,2) obey the anticommutation relation 
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where the four-momentum P M (/i = 1,2,3,4) is 
the generator of translation symmetry, which is 
broken on the latticized space-time. This led 



Curci and Veneziano |^| to propose a lattice for- 
mulation of the N = 1 supersymmetric Yang- 
Mills theory based on Wilson lattice fermions. 
The supersymmetry at zero gluino mass is then 
expected to be restored in the continuum limit, 
similarly to the axial symmetry in QCD at zero 
quark masses. The continuum limit at non-zero 
gluino mass is the corresponding theory with 
softly broken supersymmetry, which is actually 
the interesting case for phenomenology. 

The similar lattice formulation of N — 2 super- 
symmetry H is, on one hand, even simpler be- 
cause there the gluino is a Dirac particle, in con- 
trast to the N = 1 case with Majorana gluinos. 
On the other hand, the N = 2 supersymmet- 
ric Yang-Mills theory has some new generic fea- 
tures of supersymmetry, for instance, flat direc- 
tions in the effective potential and asymptotically 
free Higgs mechanism. Theoretically very ap- 
pealing is the automatically vector-like nature of 
N=2 SUSY, which is very welcome on the lattice, 
where the formulation of chiral gauge theories en- 
counters great difficulties. 

2. AN ALGORITHM FOR GLUINOS 

Although the N — 2 extended supersymme- 
try has many interesting theoretical features, it is 
better to start the non-perturbative lattice stud- 
ies with the N = 1 case, which is simpler. 

2.1. Lattice action 

The lattice action proposed in ref. ^ for the 
N = 1 supersymmetric Yang-Mills theory is 
based on the Wilson action formulated for a 
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Dirac fermion in the adjoint representation. The 
fermionic part of this action is: 

4 

a: /i— 1 

+ &v£ BM (i - 7m)^+a] } = W > ( 2 ) 

where the adjoint gauge field on links is 

K,,x, = V rs , x „[U] = 2Tr(U^T r U XIM T s ) , (3) 

and T r (r = 1, 2, 3) denotes the SU(2) generators. 
The two Majorana field components of the Dirac 
fermion field ip are defined as 



*w = -L(v, + c/), 

In terms of these we have: 



(4) 



(5) 



Therefore a Majorana fermion corresponds to 
Nf = i number of flavours. As a consequence, 
in the bosonic path integral obtained by integrat- 
ing out the Grassmann variables, the effect of a 
Majorana fermion is taken into account by the 
square root of the fermion determinant -y/det Q. 

2.2. Two-step multi-bosonic algorithm 

The numerical simulation of the N = 1 su- 
persymmetric Yang-Mills theory with dynamical 
gluinos requires the generation of -y/det Q in the 
bosonic path integral. Neglecting possible sign 
changes of det Q, we can set with Q = j 5 Q 

|Vdet(Q)| = {det(QtQ)}V4 = {det( g2 )} i/4 _ (6) 

In the framework of multi-bosonic algorithms R] , 
these can be approximated by the determinant of 
polynomials. In order to decrease the number of 
pseudofermionic boson fields, and hence the auto- 
correlation and the storage requirements, one can 
take a two-step approximation scheme The 
approximation up to a deviation norm 5 is then 



Here the polynomial P(x) is a crude approxima- 
tion of the function x- 1 / 4 in the interval x G [e, A] 
containing the spectrum of Q 2 . This approxima- 
tion is then arbitrarily improved by the product of 
two polynomials P(x)P(x). In the fermion sim- 
ulation algorithm || the second factor P(x) is 
taken into account by a "noisy correction step" , 
as introduced in ref. ||. (A similar correction 
scheme has been proposed in the multi-bosonic 
framework in ref. |10|.) For the second polyno- 
mial one does not introduce bosonic fields, there- 
fore no storage is needed. The autocorrelation is 
mainly determined by the order of the first poly- 
nomial, which can be kept low. 

The definition of the deviation norm 6 for the 
necessary polynomial approximations has to al- 
low for a flexibility in choosing functional forms, 
weight factors and approximation regions. This 
can be achieved by 



dx [1 - x a P(x)Y 



(8) 



{det Q 2 } 



2il/4 6 



det[P(Q 2 )]det[P(Q 2 



(7) 



where the function to be approximated is x~ a . 
The advantage of this choice is that S 2 is a 
quadratic form of the polynomial coefficients and 
hence the minimum can easily be found. 

Compared to Chebyshev approximations, the 
polynomials minimizing S in (||) give somewhat 
less uniform relative deviations in the interval 
[e, A]. (Remember that the Chebyshev polynomi- 
als are minimizing the maximum of the relative 
deviation, the so called "infinity norm".) How- 
ever, in the present case the greater flexibility of 
the above quadratic scheme is more important. 

3. NUMERICAL SIMULATION 

There are several interesting questions which 
can be investigated by numerical simulations in 
the N = 1 supersymmetric Yang-Mills theory. 
For a first orientation the knowledge of the gluon- 
gluino-ball spectrum of hadronic states is very 
useful. In the supersymmetric limit there have 
to be degenerate supermultiplets. For non-zero 
gluino mass these multiplets are split up, and 
for infinite gluino mass only the glueballs of pure 
gauge theory remain in the physical spectrum. 
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It is plausible (but not absolutely necessary) that 
the lowest + glueball state becomes a member of 
the lowest chiral supermultiplet, which also con- 
tains a Majorana spinor and a - pseudoscalar 
boson. This latter can be most simply built 
from two gluinos. The Majorana spinor gluino- 
glueballs are obtained as bound states of a gluino 
with gluons. 

3.1. First results on the spectrum 

The two-step multi-bosonic algorithm has sev- 
eral parameters which have to be optimized in a 
given situation (depending on bare parameters, 
lattices etc.). In ref. § a first look on autocorre- 
lations has been taken for the gauge group SU(2) 
at bare gauge coupling (3 = 2.0. (Another pos- 
sible simulation algorithm based on classical dy- 
namics type equations has recently been tested in 
ref. 0.) For an interesting spectrum calculation 
one has to go, however, to larger /3-values. 

In pure gauge theory the continuum behaviour 
starts to set in at j3 = 2.3. Therefore one can 
try first this value. For hopping parameter val- 
ues K = 0.14 — 0.15 short autocorrelations could 
be achieved on 4 3 • 8 - 8 3 • 16 lattices with the 
first polynomial P of order n — 10 — 12 and sec- 
ond polynomial P of order n = 30 — 40. The 
spectrum of Q^Q was well inside the interval 
[e,A] = [0.02,5.0]. 

Table 1 

String tension er 2 , glueball mass (M gg ) and 
gluino-glueball mass (M gg ) in lattice units. 



p 


K 




Mgg 


Mgg 


2.30 


0.140 


0.53(1) 


1.34(19) 


2.90(12) 


2.30 


0.150 


0.54(1) 


0.85(23) 


2.78(13) 



The first results have been obtained on 4 3 • 8 
lattice. The present statistics is based on 5000 
updating cycles, where a cycle consisted of a 
few standard sweeps on every field variable. To 
monitor the string tension a I determined the 
2® 2 Creutz-ratio \i = W22W11/W12 and defined 
a 2 = — logX2 — c. The results for this and the 
lowest + gluball mass (M gg ) and gluino-glueball 
mass (M gg ) are contained in table [jj The signal 
for the gluino-gluino states can be obtained from 



noisy estimators |8j , but this still needs some im- 
provement, in order to obtain M gg First estimates 
based on the study of the spectrum of Q indicate 
a critical hopping parameter K c ~ 0.165. 

3.2. Future goals 

Besides the spectrum, some other non-per- 
turbative questions can be investigated which 
have also been recently discussed in the litera- 
ture. Examples are: the phase structure based on 
exact /3- functions [Q , or the question of confine- 
ment versus screening at zero gluino mass Jl3[ . 
These problems seem to be well suited for nu- 
merical simulation studies which are feasible with 
presently available computer resources. 
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